The Pauli exclusion principle and beyond 
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The Pauli exclusion principle can be stated as inequality l^}\p\ip) < 1 for the electron density 
matrix p. Nowadays it is replaced by skew symmetry of the multi-electron wave function. The 
replacement leads to numerous additional constraints on p, which are discussed in this letter together 
with their physical implications, in particular for spin magnetic moment of a multi-electron system. 
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Recall, that the original Pauli exclusion principle 
claims that no quantum state can be occupied by more 
than one electron l|. This sets a bound on expected 
number of electrons in any given state ip 



W) < i. 



(1) 



Here p = Ylij(^\ a i a j\ i &)\' t Pi)(''Pj\ ^ s Dirac's density ma- 
trix of iV-electron system in state VP. Soon after Pauli 
discovery, Heisenberg and Dirac Q replaced it by skew 
symmetry of the multi-electron wave function 



# 6 A N H r C K 



(2) 



where the wedge stands for skew symmetric tensors built 
on the one-electron spin-orbital space Ti r of dimension r. 

In this Letter I first briefly review a recent solution 
of a longstanding problem about impact of the Dirac- 
Heisenberg replacement on the electron density matrix 
[H, Q ■ It goes far beyond the original Pauli principle and 
leads to numerous extended Pauli constraints indepen- 
dent of ([1]). In the rest of the paper I discuss some of 
their physical implications. 

The problem has a long and complicated history. A 
difference between the two versions of Pauli principle 
is best seen in Heisenberg's exchange term in electron 
energy. By adding this term Fock drastically improved 
Hartree mean field theory which initially complies the 
exclusion principle only in its original form. The ex- 
change energy plays also crucial role in formation of cova- 
lent bonds, as well as in many other two-electron effects. 
This urged Lowdin and Shull [H thoroughly analyze two- 
electron wave functions, whose structure is closely related 
to eigenvectors and eigenvalues of the corresponding den- 
sity matrix p. Lowdin Q coined for them the terms nat- 
ural orbitals and natural occupation numbers. It turns 
out that the natural orbitals of a two-electron system 

and the total state takes the 



split into pairs ip^^f 



The orbitals ^>j 1 \i/ , i nave 



2 , and 



form ^ 

the same occupation numbers A^ = A^ ' = \a% 
therefore spectrum of p is doubly degenerated. 

Beyond two electrons A 2 H r and two holes A r ~ 2 H r [7[ 
the progress was very slow. In early seventies Borland 
and Dennis Q tried out numerically the simplest such 
system A 3 H.6 and discovered the following constraints on 



the natural occupation numbers A = Spec p arranged in 
decreasing order Ai > A 2 > • • • > A r and normalized to 
the number of electrons N = Tr p: 

Ai + A 6 = A 2 + A 5 = A 3 + A 4 = 1, A 4 <A 5 + A 6 . (3) 

A rigorous proof has been found about the same time by 
M.B. Ruskai, who published it only recently 

These (in)equalities may challenge our physical intu- 
ition. They allow exactly one electron in two symmet- 
rical natural orbitals and clearly supersede the original 
Pauli principle, which in this setting reads A 4 < 1. This 
example revealed a dramatic effect of the skew symme- 
try condition on structure of the density matrix and pro- 
voked the authors' emotional comment 

"We have no apology for consideration of 
such a special case. The general 7V-repre- 
sentability problem is so difficult and yet so 
fundamental for many branches of science 
that each concrete result is useful in shedding 
light on the nature of general solution." 

Shortly afterward Peltzer and Brandstatter [l(| pub- 
lished a false solution, claiming that two-electrons, two- 
holes, and the Borland-Dennis system A 3 1i.e are excep- 
tional, and in all other cases the Pauli constraint A 4 < 1 
is the only one. After that the problem has stalled for 
more than three decades. 1 

Now, with a complete solution at hand, it may be ap- 
propriate to return back to physics and look for possible 
manifestations of the extended Pauli constraints. 

Let me first give a sample of results for a three-electron 
system A 3 7i r of even or infinite rank r: 



A fc +i+A r _ fc <l, 0<k<r; 

A 2 + A 3 + A 4 + A 5 < 2, Ai + A 3 + A 4 + A 6 < 2, 
Ai + A 2 + A 5 + A 6 < 2, Ai + A 2 + A 4 + A 7 < 2; 

Ai + A 2 + A 4 + A 7 + An + Ai 6 H < 2. 



(4) 
(5) 
(6) 



1 Mostly because of lack a proper theoretical insight, that emerged 
much latter [Til , [T2I , ll.'jl . and inconclusive results of numerical 
studies. 
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The inequalities ((4]) improve the original Pauli condition 
A < 1 . Due to normalization Tr p = 3 they turn into 
Borland-Dennis equations © for r = 6. The quadruple 
of inequalities §5§ provides a full set of constraints for 
r < 7. A proper infinite extension of the last inequality 
in the quadruple is given by ^ where the differences be- 
tween successive indices are natural numbers 1, 2, 3,4,.. .. 
Similar extensions for the other three constraints can be 
obtained by adding to them a tail of the series ((U starting 
with the term An. 

A complete set of constraints heavily depends on the 
rank r and the number of particles. Currently they are 
known for all systems of rank < 10. To give an idea of 
complexity of the problem note that for systems A 3 7iio, 
A 4 Hio, A 5 7iio the constraints amount to 93, 125, 161 
independent inequalities respectively. 

Pinned state effect. Recall, that the Pauli principle is a 
purely kinematic constraint on available states of a multi- 
electron system. Not even a minuscule violation has been 
detected so far, in spite of numerous and increasingly so- 
phisticated attempts 14| . Therefore whenever a dynam- 



ical trend is in conflict with Pauli constraints, the latter 
would prevail and the system eventually will be trapped 
in the boundary of the manifold of allowed states. This 
can manifest itself in degeneration some of the extended 
Pauli inequalities into equations. In this case the system 
and its state vector will be called pinned to the degener- 
ate Pauli inequalities. Since a pinned system is primary 
driven by Pauli kinematics, rather than by Hamiltonian 
dynamics, it should be robust and remain pinned under 
a reasonably small variation of the Hamiltonian. This 
stability gives us a chance to detect pinned states and 
in a sense makes them real. A pinned system is essen- 
tially a new physical entity with its own dynamics and 
kinematics. 

The concept can be illustrated by a string pendulum. 
For a small initial horizontal impulse at the lowest point 
its mass is confined to a circle. However, if the energy is 
big enough to raise the mass above the suspension point 
but is not sufficient to neutralize gravitation at the apex, 
then at some moment the string constraint becomes ir- 
relevant and the mass switches into a free trajectory. 

As a more relevant example, consider the first excited 
state of beryllium atom of spin S = 1 and S z = 1. Here 
are its natural occupation numbers calculated from 10 
spin-orbitals in the lowest three shells Is, 2s, 2p [IB] ] 



1.000000, 0.999995, 0.999287, 0.999284, 0.000711, 
0.000707, 0.000009, 0.000007, 0.000000, 0.000000. 



(V) 



Observe that the first orbital is completely filled and the 
last two are empty, i.e. the state is pinned to the initial 
Pauli constraints < A < 1. This kind of degeneration is 
very common. In its ultimate form it has been elevated to 
the Aufbau principle and incorporated into Hund's rules. 

The completely filled and empty orbitals are inactive, 
and therefore we effectively deal with a reduced system 



/\ 3 H.7, in which all orbitals are partially filled. The ex- 
tended Pauli constraints for this system amount to the 
quadruple of inequalities ([5|). Note that the exclusion 
principle itself Ai < 1 doesn't enter into the list explic- 
itly, but can be derived by taking sum of the second and 
third inequalities that gives Ai + Xq — A7 < 1. One can 
infer from this that the Pauli degeneration Ai = 1 for 
a system of this format would pin it to the last three 
constraints ([5|), that still fall short of Ai = 1. Actual cal- 
culation for beryllium data ([7|) shows that one of these 
three inequalities turns into exact equation 



Ai + A2 + A4 + A7 



2, 



(8) 



and the other two are exact within a rounding error. 

This equality stretches the Pauli principle to its limit, 
and therefore restricts the system's kinematics. Recast- 
ing it into the form 



(ajar 



a\a 2 



a\a4 



a\a 7 )^ = 2* 



(9) 



we get a selection rule for Slater determinants [i,j,k] = 
if>i A ipj A ipk that enter into the decomposition ^ — 
^2 Cijk[i, j, k]. Such a determinant should include two 
natural orbitals from the set {1,2,4,7}, and one from 
its complement {3,5,6}. As a result the pinned system 
splits into two components A 2 H.4 (8 H3, spanned by the 
above quadruple and triplet, and containing respectively 
two and one electron. The components can be readily 
identified with two spin-up electrons, and one spin-down. 

The selection rules coming from the other two pinned 
inequalities bound the decomposition to four Slater de- 
terminants 

= a[l, 2, 3] + /3[1, 4, 5] + 7 [1, 6, 7] + S[2, 4, 6] (10) 
with the amplitudes, up to phase factors, set by equations 
' 2 = A 3 , \(3\ 2 = \ 5 , M 2 = A 7 , 



W = *3> 
M 2 + H 2 = A 2 , 



|5| 2 = A 4 , | 7 | 2 + H 2 = A 6 . 



This gives us an insight into structure of the pinned beryl- 
lium state whose form (fT0| should be stable under a per- 
turbation. It corresponds to the type VII in Schouten 
classification of 3- vectors \I/ € A 3 H.7 16j |. 

A similar consideration can be applied to the Borland- 
Dennis system A 3 7ie in which every state is pinned by 
two equations ([3]). They reduce its internal kinematics 
to that of three qubits H2 <8> H2, each spanned by 

a couple of symmetric orbitals tpi^j, i + j = 7. In this 
interpretation the respective pairs of occupation num- 
bers Xi,Xj become spectra of three reduced states of 

€ H2 ® H2 ® Ti-2, and the inequality ([3]) turns into 
Higuchi-Sudbery-Szulc compatibility condition [17| . 

Spin-orbit interaction versus Pauli constraints. In the 
above setting we were unable to address specific spin ef- 
fects buried in join spin-orbital space Ti. r = Jig ® TL S . 
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Actually the orbital Hi and spin H s degrees of freedom 
play very different role. The former, via Coulomb in- 
teraction, are primary responsible for dynamics, whereas 
the latter, disregarding a small relativistic correction, are 
purely kinematic. In this approximation the total spin is 
a good quantum number, that single out a specific com- 
ponent in Weyl's decomposition 



(He ® Ti s ) = Ke^K* 



(11) 



\u\=N 



Here H\ is an irreducible representation of the orbital 
group SXJ(Hi) defined by Young diagram v and H v s is 
the representation of spin group SU(W S ) with transpose 
Young diagram v*. The electron spin group is SU(2) 
and in this case v is a two-column Young diagram. For 
example, in the ground closed-shell state of beryllium 
atom v — , while for the excited state of spin one 

v = gP. In general the total spin S is equal to half- 
difference of lengths the two columns of v. 

For spin resolved state \& £ H^ <g> H v s the Pauli con- 
straints amount to inequalities between orbital \ and 
spin /ij natural occupation numbers normalized to traces 
iV" and 1 respectively Q . As a toy example consider three 
electrons in c£-shcll (dim 7-^ = 5) in low spin configuration 
v = pp where the constraints are as follows j4j j , 6. 1 ] 

Ai + |(A 4 + A 5 ) <2, 

H < 3 - 2(Ai - A 2 ), m<3-2(A 2 -A 3 ), 

fi > 2(Ai - A 3 ) - 3, n > 4Ai - 2A 2 + 2A 4 - 7. 

Here ji = H\ — /i 2 is spin magnetic moment in Bohr mag- 
netons hb and I set g = 2 for the electron gyromagnetic 
factor. 

The first inequality in the list improves the Pauli con- 
dition Ai < 2 for the number of electrons in an orbital. 

The second line provides a kinematic upper bound on 
the spin magnetic moment. This contrasts sharply with 
prevailing opinion that the reduction of magnetic mo- 
ment is a dynamical effect caused by spin-orbit coupling. 2 

One has to be careful with the lower bound for ji = 
\i\ — /x 2 given in the last line. Recall, that fJ-i,fJ-2 are 
probabilities to find the system in spin-up and spin-down 
states with respect to some natural direction determined 
by leading eigenvector of the spin density matrix. It may 
not coincide with the direction of magnetic field H used 
for measurement of the moment. However, for negligible 
spin-orbit coupling the natural direction is free to orient 
itself along the field to produce the maximal moment 
fi = fii — H2- Otherwise the coupling might reduce this 



ideal value to the level that breaks the Pauli constraint. 
This would be a signature of an appreciable spin-orbit 
interaction. 

This example provides a key to an old puzzle about 
reduction of magnetic moment of a transition metal atom 
in a crystal, relative to its value in free space. I will focus 
below on body-centered cubic form of iron (a-Fe), stable 
at normal condition, using a recent high precision 7-ray 
diffraction study by Jauch and Reehuis 1131. 

The study support an earlier finding [20|, bases on 
measurement of a Mossbauer shift, that the electron con- 
figuration of 3<i-shell in iron crystal is d 7 rather than d 6 as 
for free atom. Such a modification of metallic 3d-shell is 
pretty common, albeit not universal. It reduces the max- 
imal spin magnetic moment per atom from 4 /ib to 3 \ib-, 
while the experimental saturation moment is 2.22(1) \ib- 

Recall, that in a crystal field of cubic symmetry the 
orbital moment is quenched and d-shell splits into two ir- 
reducible components of dimension 3 and 2 called i 2s and 
e g subshells. The orbital density matrix retains the crys- 
tal symmetry, and therefore reduces to scalars nt and n e 
on the above subshels. As a result the orbital occupation 
numbers A = (nt,nt,nt,n e ,n e ), 3n* + 2n e — 7 depend 
only on one parameter nt, where I expect n t > n e . 

The Pauli constraints on spin magnetic moment /u, ver- 
sus the orbital occupation number nt for high spin d 
configuration are shown in Fig. [TJ The plot also includes 
the experimental data for iron /1 = 2.22, n t — 1.458(7). 
The latter value is deduced from the relative number of 
d-electrons 62.5(3)% in i 2ff -subshell 

The data suggest that the observed magnetic moment 
is the maximal one kinematically allowed by the electron 
density, i.e. the iron (i-shell is pinned to the boundary 
segment [A, B] given by equation /i = 7n t — 8. Let me 
stress again that the constraint on spin magnetic moment 
M ^ 7rit — 8 comes from the kinematics of e?-shell and has 
nothing to do with spin-orbit coupling. One can interpret 
it as an ultimate result of a multi-electron "exchange in- 
teraction" , though the term suggests a dynamical origin 
of the effect, while in fact it is purely kinematic. 

The constraints shown in Fig. [1] are actually derived 
from a system of 55 inequalities between orbital A^ and 
spin Hj natural occupation numbers. The pentagon 
ABCDE is a projection of a multidimensional polytope 
that includes all spin occupation numbers fa, rather than 
just the total magnetic moment fi = 3/^i + /z 2 — ^3 — 3/X4. 
It turns out, however, that the boundary segment [A, B] 
has unique pull back [A, B] in the bigger polytope 
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2 See for example fTsL p. 98]: "If spin-orbit interaction is zero free 
spin magnetism is realized. This is meant by saying that the 
orbital moment is quenched." 



where the first 5 components represent the orbital occu- 
pation numbers and the rest are spin ones. This allows 
to recover spin occupation numbers fa of iron d-shell 



(0.69, 0.23, 0.08, 0). 



(12) 
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The small probability 0.08 to have negative magnetic mo- 
ment — 1 /ib can be seen in a spin density plot [22], p. 170] 
built from polarized neutron scattering data. I wonder if 
the data contain enough information to recover the spin 
occupation numbers and verify the values (|12[) . 

Magnetovolume effect. Pressure is an effective tool to 
change electronic structure, though apparently it remains 
detached from spin. This makes the observed reduction 
of magnetic moment under pressure a mystery, known as 
magnetovolume effect 2l|. I use it below to test robust- 



ness of pinned iron states. 

For a-iron the magnetic moment is governed by the 
pinning equation [i — 7n t — 8, at least when the pres- 
sure is not too high. If at some point the constraint 
would cease to exists, then the frozen degrees of freedom 
will be released and change the dynamics of the system. 
No such change is observed until a transition from bec 
to hep crystal structure has started. The magnetic mo- 
ment at this point is about 1.9 [Ib [H[ and the electron 
density should be close to spherical one n t ~ n e ~ IA 
to accommodate the two incompatible symmetries. The 
transition actually passes through a phase coexistence 
interval. A linear extrapolation to the point where bec 
phase vanishes gives the magnetic moment about 1.8 [Ib- 
This is consistent with the evolution along the pinning 
edge [B,A]. 

The effect is more tricky in a face-centered cubic form 
of iron stable at high temperature or low density. It has 
two competing spin states, first conjectured by Weiss [23| . 
The one with high moment /x 2.5 \ib at some critical 
pressure collapses into a low spin one ftsil \xb ■ 

Recall that 2.5 /is is the maximal magnetic moment 
kinematically allowed for cubic iron, provided n t > n e . 
Assuming the latter condition we infer that the high 
spin state should be pinned to the point B and remains 
trapped at this point until one of the two pinning equa- 
tions /i = 7n t — 8 or /i = 16 — 9n t is released. As the above 
scenario suggests, at some critical pressure the first con- 
straint is eventually freed and the state evolves along the 
segment [B, C] into the low spin state C, where it will be 
stabilized until the second pinning constraint /i = 16— 9nt 
is relaxed. Then it slides along the path [C, D] into a non- 
magnetic phase. This picture qualitatively fits the cal- 
culated 24 1 and the observed [25j evolutions and gives a 
hint about possible origin of the two spin states. 

May be I have to conclude with a disclaimer that this 
letter is not about iron magnetism, used here as a testing 
ground for new concepts and methods originated from a 
recent advance in understanding of quantum kinematic 
constraints. Until now the kinematic effects were largely 
neglected or treated as dynamical ones. The confusion 
was clearly recognized by chemists [26j], and the above 
examples may support their objection. In the end, the 
role of Hamiltonian in some cases may reduce to pinning 
a state to a certain kinematic constraints which qualita- 



B=(3/2, 5/2) 



A=(7/5, 9/5) 



C=(5/3, 1 ) 



E=(7/5, 0) 



FIG. 1: Pauli constraints on spin magnetic moment (ub) for 7 elec- 
trons in ti-shell in cubic crystal field versus the occupation number 
nt of a t2 g orbital. All points within the pentagon ABODE are 
admissible. A black dot represents experimental data for iron. 



tively shape the system's behavior. 
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